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Abstract 

For basic discrete probability distributions, — Bernoulli, Pascal, Poisson, hypergeo- 
metric, contagious, and uniform, — g-analogs are proposed. 

1 Introduction 

g-analogs of classical formulae go back to Euler, g-binomial coefficients were defined by 
Gauss, and g-hypergeometric series were found by E. Heine in 1846. The g-analysis was 
developed by F. Jackson at the beginning of the 20*'* century, and the modern point of 
view subsumes most of the old developments into the subjects of Quantum Groups and 
Combinational Enumeration. 

The general philosophy of q-analogs is that of a deformation, with the deformation 
parameter q being thought of as close to 1. This point of view is certainly not all- 
encompassing; for example, representations of Quantum groups when q is a root of unity 
are of independent interests; more importantly, the g-pictures sometimes possess proper- 
ties singular in {q — 1) or otherwise not regularly dependent on (g — 1); regularization of 
divergent /infinite {q = 1)— quantities is another useful feature of g-analogs... the list goes 
on. 

The typical example is 

lim [x] = , \q\ < 1, (1.1) 



X— >oo 



where 



1 



(1.2) 



is the g-analog of a number (or object) x. (A quick introduction to the g-calculus is 
available from many sources, e.g. Chapter 2 in [6].) 

More examples of such sort will be found below in this paper, the 1'^* one in a series 
devoted to a ^-probabilities. In the next 6 sections we look at g-analogs of Bernoulli, 
Pascal, Poisson, hypergeometric, contagious, and uniform distributions, respectively. It is 
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surprising how many new effects appear compared to the classical theory at g = 1. For 
example, even for the Bernoulli distribution, the profound classical differences between 
finite and infinite number of trials are mitigated when q enters the picture, so that one 
can write down the probability (formerly zero) of many individual events of infinite type, 
such as 



{i-pr= n~o(i-pg'), |9|<1, 



(1.3) 



the probability of coming up with all "tails" during an infinite number of coin flips; the 
probability of coming all tails during n coin flips is 



(1.4) 



More generally, the probability of coming up with precisely k heads during an infinite 
number of coin flips is 



(1 -(/)...(! -g- 



■(l-p)°°, kGN. 



(1.5) 



The six probability distributions discussed in this paper are all discrete, univariate, 
basic, and relatively simple. More classical distributions can be found in [1, 3-5]. 



2 q-binomial distributions 

Suppose we have a random variable C, — "2-sided coin" — which takes two values: 1 with 
probability p, and with probability 



p' = l-p. 

After n throws, the total sum accumulated, 

= Cl + ••• + Qni 

obeys the Bernoulli distribution 

Pr{in = k) 

PiS a g-analog of this distribution we take (with < g < 1) 



n 



pi^pin-n^ < K < n. 



Pr{in = M) 



p^il -p)"-'', 0<K<n, 



where 



[x\...[x — 

= 0, s G Z+, 



, K e N, 



(2.1) 
(2.2) 

(2.3) 

(2.4) 

(2.5a) 
(2.5b) 
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are the g-binomial coefficients, and 

[0]! = 1; [k]\ = [1]...[k], «gN, 

are the g-factorials. 

To justify formula (2.4), we need to prove that 



(2.6) 



K 

k=0 '- 



n 



(2.7) 



This formula follows from the following identity: 



K=0 



K=Q 

when V = p; here 

(a+6)" := n^-^ (a + q'b), n G N, {a+bf := 1. 
Formula (2.8) is, in turn, the b = 1-case of the general formula 



Y: I a-ib+vr- = E 



K=0 



K 
K=Q 



n 



(2.8) 



(2.9) 



(2.10) 



b^ia+v)'^-". 

To prove formula (2.10), let us use the easily checked by induction on m Euler's formula 

(2.11) 



{x+yr = E 

3=0 



m 
IJ J 



Then the LHS of (2.10) can be rewritten as 



n 




n — K 






i 







(2.12L) 



while the RHS of (2.10) can be similarly rewritten as 



E 

s,e 



n 




n — s 




b' 


s 




£ 



a— Vg(2), 



(2.12R) 



and these two double sums bijectively coincide, for each fixed i, when s is identified with 
n — £ — K. 

Now, to calculate the expectation values of powers of we use the easily proved by 
induction on m G N formula 



X 



_d_ 

dqX , 



)"^|:i;^(|:[^:']Mr.«i«-»r-').<^) 



(2.13) 
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where 



df f{qx) - fix) 



dqX 



{q — l)x 
is the g-derivative: 
d 



dqX 



In particular, for m = 2, we get 



X 



dqX 



X- h qx 

dqX 



dqX 



Applying the operator ( p- — 

V dqp 

n 

< 6^ >= E{Q = 

k=0 

n 



s = 1, 2, to formula (2.8), we obtain 



v=p 



K = 



p-{l -pT-- = [n]p, 
p'^{l-pr-'' = [n]p + qp^[n][n-l], 



where we used the obvious relation 

Prom formulae (2.17) and (2.18) we find that 

VariU) =<?n>-< Cn = [n]p{l - p). 



(2.14) 



(2.15) 



(2.16) 



(2.17) 
(2.18) 



(2.19) 



(2.20) 



Notice that formulae (2.4), (2.17), (2.18), (2.20) have a well-defined limit as n ^ oo: 



Pri^oo = M) 



P 



[i^]\\l-q 



(1 -pr, 



< Coo >-- 



^ Soo ^ 



P 



1-q 
P 



1 

1-q 

^ P(l - P) 
1-q 



P 



1-q 



(2.21) 
(2.22) 

(2.23) 
(2.24) 



So far, we have treated the random variable ^„ as an object in its own right. Let us now 
turn to the representation of as a sum of n "coin" throws, as reflected in the classical 
formula (2.2). This sum-formula (2.2) remains intact under g-deformation. However, for 
general q, the random variables Q^s are no longer independent or identically distributed. 
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1, with probability p, 25) 
0, with probability 1 — p. 

For K G N, suppose the random variables (i, ...,Ct have already been defined. Denote by 

<^K = ctniCii ■■■iCk) the number of zeroes appearing among the values of the random vari- 
ables (i, The random variable Ck+i takes the values 0, g**, with the conditional 
probabilities 

Pr(CK+i = 0|a« = r) = 1 - q'p, < r < k, (2.26a) 
Pr(a+i = q%K = r) = Si_,q'p, < r < k. (2.26b) 

For example, 

^^^((2 = 0|Ci = 0) = 1- qp, (2.27a) 

Pr{C2 = 0\Ci = l) = l-p, (2.27b) 

Pr{C2 = l|Ci = 0) = qp, (2.27c) 

^r(C2 = g|Ci = l)=P- (2.27d) 

By induction on k, it is easily seen that if the last before non-zero value appearing 
among (i, ...(^^ was q^, then Ck+i can take only the values and q^^^ with non-zero prob- 
abilities; if all the Ci5---;Ck took value 0, then Ck+i takes only the values and 1 with 
non-zero probabilities. 

A better description of the same distribution is possible, if instead of conditional prob- 
abilities we work with joint ones. Denote by C the event of r in a row appearances of 
zeroes, r G Z+, and similarly by 

Qa(0)^OQa(l)^l __Qa(fe-l)^«-lQa(K)^ a(-) G Z+, (2.28) 

the event of a(0) of zeroes followed by g'^ = 1 followed by a(l) zeroes ... Now set 

Pr(0«(o)gO __^K-iQa(«)^ := (1 -p)2S«W n'^-^ (pgSo«W); (2.29) 

for K = 0, formula (2.29) is to be understood as 

Pr(0'^) = (1 -pf, a EN. (2.29') 

Let us now verify that the "microscopic" formulae (2.29) imply the "macroscopic" 
formula (2.4). Denote 

i 

Hi)\ = J2a{s). (2.30) 

s=0 

We have to verify that 



{l-p)''-^p^ n^^Jo^gl^Wl 

\a(K)\=n—K 



p^{l -p^-", (2.31) 
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which is equivalent to the g-counting formula 



\a{K)\=n—K 



(2.32) 



(For g = 1, we recover the classical result: the number of solutions in nonnegative integers 



n 



of the equation a(0) + ... + a{K) = n — k, is ^ J. ) 

We shall prove formula (2.32) by the double induction on N := n — k and k, in the 
form 



|a(K)|=7V 



N + K 
K 



Now, 

n«-Qlgl««l = g«(0)g«(0)+a(l)| __^a(0)+...+a(K-l) ^ qEoik-s}ais) 

Therefore, the identity (2.33) becomes 



J2 q-^o^<') =q-^k 
|a(«)|=Ar 



N + K 
K 



(2.33) 



(2.34a) 
(2.34b) 



(2.35) 



For N = 0, formula (2.35) becomes 1 = 1 no matter what k is. For k = 0, formula (2.35) 
is true by definition (2.29'); for k = 1, formula (2.35) becomes 



a{0)+a{l)=N 



-a(l) ^ ^-N 



N + 1 
1 



which is obviously true for all A''. Suppose formula (2.35) is true for the pairs (k, N = n) 
and {k — 1, N = n + 1). Consider the pair {k, N = n + 1). Let's divide the set of the a's 
with \a{K)\ = n + 1 into two groups: those with a(«;) > and those with a{K) = 0. For the 
1^* group, the set a(0) = a(0), a{K — 1) = a{K — l),a{K) = a{K) — 1, satisfies |a(«;)| = n, 
so that, by the induction assumption. 



q- ES = q-ES ^»(^)-« = q-'^q- 

\d{K)\=n 



n + K 

K 



(2.36) 



The 2"^^ group has effectively k — 1 a's, so that again, by the induction assumption. 



^g-ES"' Ws) = ^-(n+l)(«-l) 

We thus have to check that 



n + K 

K-1 



-K(n+1) 



n + K 



(n+l)(K-l) 


n + K 


= g-(n+l)K 


n + 1 + K 




K-1 




K 



(2.37) 



(2.38) 
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which is equivalent to 



n + K 

K 



+ 



n + K 

K-1 



n + 1 + K 

K 



(2.39) 



which is obviously true. 

Notice that for a^n) = oo, formulae (2.29) become 

Pr(0"(°)g°...g'^-^0°°) =p''{l -p)~gEo"'(«-^)«W, 
Pr(0°°) = (1 -p)°°. 



(2.40) 
(2.40') 



As in the classical theory (cf [10] p. 59), we can calculate the probability of observing 
< K {< n) zeroes in n trials: 



Pr{an < k) = ^ 



1=0 



n — 1 

K 



Jo 



p--(l -pY = [n] 

= ( r x'^-^-'^{l -qx)^dgx)/{ /^x"-^-''(l -qx)''d, 
Jo Jo 

Similarly, the probability of < £ (< n) non-zeroes in n trials is 



X). 



(2.41) 



Pr{an > n - ^) = 



i=0 



p\l -pT-^ = 1 - [n] 



n — 1 
£ 



r x^l-qxf-^-^dqX. 
. Jo 



(2.42) 



(Formulae (2.41) and (2.42) can be easily proven upon (/-differentiation with respect to p 
and using the obvious relation 



d 

dqp 



(l-p)" = -[n](l-gp)"-^ 



(2.43) 



In the limit n — > oo, this probability becomes 



Many, if not all, classical formulae in probability have g-analogs. Let's take a look at a 
few of such formulae involving higher moments for the Bernoulli distribution. 

d \^ 

to formula (2.8) and using the relation 

v=p 



First, applying the operator p^ — — 

\dqP 



X = K =^ K — I 



r{x - [z]), 

we get 

E{nt^q-\in - m = P'- n[-o^ [n - i], 
a g-analog of the familiar formula 

< Cn(Cn - l)-(Cn - r + 1) >= /n...(n - r + 1). 



(2.45) 
(2.46) 
(2.47) 
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Second, let 

/x;=<C>, Mr =< (Cn- < >r >, ^ G Z+, (2.48) 

be the moments, around zero and < > respectively, considered as functions of n, r. 
Romanovsky [9] has proved that 

/x;+i = {np + p{l-p)-^M), (2.49) 

fJ'r+i=Pi'^-p){nriJ,r-i + ^)- (2.50) 

Formula (2.49) has the following g-analog: 

= {[n]p + p{l-p)^){4). (2.51) 

Formula (2.50) has no clear g-analog, in part because the notion of higher central moments 
is not unique in g-probability. Certainly, the classical definition 

fir = E{{C- < C >)'■) (2.52) 

is not useful, as the objects 

-^iia + PxY) (2.53) 

dqX 

lie outside compact formulae of q'-analysis. Two other possible definitions are 

fXr = E{{C-<C >r), (2.54) 



where [7] 



n 



K 



a'^b''-'', n e Z+, (2.55) 



«;=0 

and 



fir{s)=E{{^-q^<C>Y). (2.56) 
Using the definition (2.56) for the g-Bernoulli distribution (2.4), we find 

Hr+i{0;p) = p{l - p){q[n][r]i^r-iihPQ) + -j-il^rihp))), (2.57) 

dqP 

Hr+i{-r;p) = p{l - p){q^''[n][r]fir-ii-r;pq) + — (/Xr(-r;p))). (2.58) 

These formulae indicate that the central limit theorem for the g-Bernoulli distribution 
may not exist, at least in the classical sense. 
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3 q-analogs of negative binomial distributions 



The negative binomial distribution, also called Pascal distribution, can be arrived at via 
many different routes. Perhaps the simplest one is as the waiting time in a succession of 
Bernoulli trials until the appearance of r*'* non-zero for the first time. 

Let's first consider the case r = 1. Let W be the random variable, waiting time until 
first non-zero. By formula (2.29), 



Pr{W = j) = Pr{0^-\°) = (1 -py~\^~^p, j G N. 
We further g-modify this formula by setting 

Pr{W = [j]) = (1 -py-'q'-'p, j e N. 
This is our q'-analog of the geometric distribution. Since 

oo 

^(1 -py~\^-'p = 1 - Prion = 1 - (1 -pr, 

we get a g-analog of the formula for the sum of a geometric progression: 
1 - (1 -p)°° 



(3.1) 
(3.2) 

(3.3) 



s=0 



P 



(3.4) 



Remark 3.5. Most of the formulae appearing in this paper remain true when q is con- 
sidered as a formal variable, or as a complex one (with occasional restrictions of the type 
\q\ < 1, \q\ > 1, etc.) It is only for the sake of probability interpretations that q is 
considered to be a real number between and 1. 
Formula (3.4) has a finite counterpart. 



N 



X) q 



(1-x) 



N+l 



(3.5) 



s=0 



This relation is easily checked by induction on N. 

For general r G N, the probability that the r*'* non-zero occurs at exactly the j^^ trial, 
j > r, is, by formulae (2.29), (2.30), (2.33): 

Pr(0"''g°...0"('-^)g'^-^) = ^(1 -p)J-Vg''^^"'"^g"^o~'^"^'^ 

\a(r-l)\=j-r 

j — r + r — 1 



(1 ^py-^^pr qr{j-r) ^-{j-r){r-l) 



r-1 



(1 -py-w-'-p' 



r-1 



(3.6) 



Since the probability of having exactly i non-zeroes during an infinite number of trials is, 
by formula (2.21), 



1 



P 



(1 ^P)°°, 



(3.7) 
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the probability of having < r non-zeroes is, therefore, 



''""^ 1 / p 



(1 -pY 



(3.8) 



Thus, 



E(i -py 

s=0 



r — 1 + s 




s 





(3.9) 



(- 



r-1 



This identity can be gotten directly from formula (3.4) by applying the operator 
to it. 

Notice that formulae (3.3) and (3.9) show that our g-distributions do not sum up to 1 
and therefore have to be re-scaled. For example, formula (3.2) becomes 



(3.10) 



Remark 3.11. Formula (3.6) could have been arrived at via one of the standard routes 
as the conditional probability of having r — 1 non-zeroes at the I''* j — I throws, with 



probability 



1 



p^ (1 —pY ^ by formula (2.4), followed by a non-zero appearing at 



the j^^ throw, with probability q^~^p by formula (2.26b). 

We conclude this section by re- visiting the "probleme de parties", one of the first 
problems in probability discussed and solved by Fermat and Pascal in their correspondence. 
In essence, we want to find the probability that a non-zeroes appear before b zeroes in the 
Bernoulli trial of a + 6 — 1 throws. This can happen in either one of the b ways, when the 
a*'* non-zero appears at the {a + £Y^ trial, < ^ < 6 — 1. By formula (3.6), the probability 
of this is {j = a + £, r = a): 



(1 - pYq'p- 



a + l-1 
a-l 



so that the total probability is 

b-i 



p 



a + i-l 
a-l 



(1 -pYd'. 



(3.12) 



(3.13) 



Similarly, the event that b zeroes appear before a non-zeroes can happen in one of the a 
ways, when the 6*'* zero appears at the (6 -|- £)*'* trial, < ^ < a — 1. The probability of 
this event is, by formula (2.26a): 



b + £-l 
6-1 



p\l-pY-\l-q'-'p) 



6 + £- 1 
6-1 



p%l -pY. 



Thus, the total probability of b zeroes appearing before a non-zeroes is 

a-l 



P2 = {1-pYY1 



b + e-l 
6-1 



P 



(3.14) 



(3.15) 
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Since Pi + P2 = 1, we find 

6-1 r 



=1 '- 



a + i-1 
a-1 



a-l 

(i-p)V + (i-p)'E 



b + i-l 
6-1 



/ = 1, Va, 6 G N, 



(3.16) 



an identity which is not immediately obvious. 

The same probability Pi can be calculated differently, as the outcome, out of a + 6 — 1 
trials, of a + s non-zeroes, < s < 6 — 1. Thus, 



s=0 



a + b-1 
a + s 



p^+'il -p) 



b-l-s 



(3.17) 



and wc arrive at another nonobvioTis (even for q = 1) identity 



a + e-1 



(1 -pYq^ = y: 



s=0 



a + b-1 
a + s 



p"+^(l -p) 



b-l-s 



(3.18) 



4 q-Poisson distribution 

One of the shortest derivations of the Poisson distribution consists of considering, as 
Poisson originally did, the limit 



n 00, pn ^ X 
in the Bernoulli distribution: 



\kJ Ki n** n 



(4.1) 



(4.2) 



The g-picture is more interesting. First, for \q\ < 1, the expression Pr{^n = M) (2-4) has 
the ra — > 00 - limit (2.21): 



1 



P 



(1 -pY 



[k]\yi-q^ 

We can get a g-Poisson distribution from this by setting 
A 



p = lim 

so that 

Pr{X = [k]) 

and therefore 



A(1-(Z), 



{1-X{1 -q)r, 



EoiX)=E,iX;q)=^^^ = ^^ 



(4.3) 



(4.4) 



(4.5) 



(4.6) 
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the well-known formula. Here 



K=0 ^ ' 

is the ^'-family of exponentials: 



dqX 

By Euler's formula (2.11), 



j=0 

so that 

Eo{X)Ei{-X) = 1; 



(4.7) 



(4.8) 



(4.9) 



(4.10) 



the latter 2 formulae are of course classic. 

Formula (4.10) can be generalized, as follows. Consider the probability generating 
function for the g-Bernoulli distribution: 



k=0 

n 

K 

n 

K 



K=0 



^ I {zpY{1 -pT-'^ [by (2.8)] = ^ 



{zp -pY 



p'^iz -Vf. 



As n — > oo, this generating function becomes 

On the other hand, since 

(1 -a;)"+^ = (1 -xY{\ -q'^xf, 
we have: 



(1 -x) 



{l^q-Pxy 



(4.11) 



(4.12) 



(4.13a) 



(4.13b) 



and therefore 

(1 -pY-'^ = (1 -pYil -q^pY'" 

{l-p)°° = Ei 



(1 ^pY 



(1 — '^pY 
P \ 1 



1 - q) EoiX) ' 



g-Probability: I. Basic Discrete Distributions 



85 



so that 



fb.,oo = lim -pT'^ = E (t^YA(^ -py 

n->oo-^^ \_k} ^-^\l — qj [k\\ 

Thus, 

which can be equivalently rewritten as 

Taking the hmit n — > oo in formulae (2.17), (2.20), we find: 

<X>=X, 

Var{X) = A(l - (1 - q)X). 



EojXz) 
Eo{X) 



(4.14) 



(4.15) 



(4.16) 



(4.17) 
(4.18) 



Let us now consider the case when \q\ > 1, so that we are walking outside the traditional 
probability theories. Again, set 



P 



X 



n 



Then, as n — > oo. 



^ [n]...[n-K + l] A- ^ A;^ (.) 



since, as n ^ oo, 

[n-e] _ g"-^ - 1 _f, 
[n] ~ g" - 1 ' ^ 

Next, by formula (4.13b), 

(1 -p)-« = (1 -T^)"-'^ 



(1 • ^ 



(1 -5«-«^)« ■ n-ci>'^ [n] 



lim (1 -X)-". 



(4.19) 



(4.20) 



(4.21) 



(4.22) 



Thus, 

Pr{X = [k]) 
In particular, 

lim (1 - A) 

n-*oo n 



-© A 

-^A''(l -A)-« lim (1 - — )", \q\ > 1. 



«:=0 



(4.23) 



(4.24) 



this formula is not a (/-analog of anything classical. 
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5 q-hypergeometric distribution 

Imagine that we have an urn consisting of two types of balls: m marked '1' and u marked 
'0'. We pick out at random one ball, record its value and leave it outside the urn; then 
proceed again, for a total of n draws. Had we returned each picked ball back into the 
urn, we would have the Bernoulli trials; since we don't return the balls, we get something 
different, called the hypergeometric distribution: the probability of ending up with k '1' 
balls out of n draws is 



u 

n — K 



N = m + u, 



(5.1) 



see [4, 5]. 

As a gr-analogue of this distribution we set 



To justify this definition we have to verify that 



m 




u 


g("^-«)(" 


K 




n — K 





N 
n 



m 




u 






_n — K_ 



Jm—K) (n— k) 



m + u 
n 



This identity results by picking the a;"-coefficient in formula (4.13a): 

(1 -2;)™(1 -q'^xT = (1 -x)™+", 
and using the Euler formula (2.11): 



(5.2) 



(5.3) 



(5.4) 



n k+£=n 



m 




u 


K 







m 




u 






_n — K_ 



{m—K) (n— k) 



m + u 
n 



where we used the obvious relation 



2 



+ 



(-a;)^©, (5.5) 



(5.6) 



Similar to the Bernoulli case, we can treat the g-hypergeometric distribution (5.2) as 
a macroscopic object and inquire about its microscopic representation. The latter can be 
guessed from the relations 



Pi^n = [n]) 
Pri^n = 0) -- 



m 


1 


'N' 


n 




n 



N 
n 



[m] [m — 1] [m — n + 1] 
[n] [A^-l]"'[iV-n + l]' 
[u] „m - 1] „m [u-n + 1] 



[A^]^ [A^-1] ■■■[A-n + 1] 



(5.7a) 
(5.7b) 
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which suggest that in the representation 

en = Cl + -+Cn 

of n successive draws, we should set 

= ("^-^)»(^) n^-o' [m - i] n]i^-' [u - j]/ n:;zl [at - 7]. 

Pr(0") = q 



m 

K 



u 

n — K 



N' 
n 



u 


1 


'N' 


a 




a _ 



(In the C-language, we have 

Pr(CVO) = M, Pr(C- = 0) = ,-M 



(5.8) 

(5.9a) 
(5.9b) 
(5.9c) 

(5.10) 



at each ball pick-out when the number of marked 'nonzero' balls in the urn is m, and the 
number of those marked 'zero' is u = — m.) 

To prove that microscopic formula (5.9a) implies the macroscopic formula (5.2), we 
need to verify that 



\a(K)\=n—K 



(5.11) 



and this equality follows from the already proven formula (2.35): 



E 



7E 5 



E' 



r™E «(s)o-E = m{n-K,) -(n-K)K 



\a{K)\=n—K 



n — K + K 

K 



(m— K)(n— k) 



In the classical case g = 1, we can rewrite formula (5.1) as 



Pri^n = k) = 
so that in the limit 



n 



n 



fe-i 



m 



n 



n—K+l 



u-e 



(5.12) 



(5.13) 



formula (5.12) becomes the Bernoulli one; the classical explanation is that when m and u 
are large, it makes little difference whether the picked-out balls are returned back to the 
urn or not. 

For g 7^ 1, the situation is more interesting. Certainly formulae (5.13) are not the 
correct ones. We proceed as follows. 
Let |g| > 1. Set 



AT ^ 00, 



[m\ 



p, 



(5.14) 



and re-write formula (5.2) in the form 
Pri^n = M) = 

Now, 



1 / [m-s\q- 
"^=0 { [iV-(n-At)-s].^ ^=0 



[N-i] 



[m — s]q 



k—n 



[-s] +q "[m] 



[N -{n- k)-s] [K-s-n] + g«-^-"[Ar] 



[by (5.14)] 



K~n_Q ^ 



\rn\ 



\u - t[q"' _ [N - m - (]q"' _ [N - (] - [m] _ _ _ 

[N-e] ~ [N-£] ~ [N-e] ~ ~l-£] + 



so that 

Pri^n = M) 



as desired. 

Suppose now that |g| < 1. Denote that q by Q. Set 
q = Q-\ \q\>l. 

Since 

[k]\Q = q<D[n]\,, 



M-a) 



Q 



we have 

Pr{in = Mq) = Q('"-«)("-'^) 

where 

k' = n — K. 
We see that our formula (5.2) in the form 



K){n—K) 


m 




u 


J 


'N' 




K 


Q 


n — K 




. ^ . 



u 




m 




'N' 






_n — k' 


J 


. ^ . 



Pr{U = k) 
allows the symmetry 



m 




u 






Jl — K_ 



(m— K)(n— k) 



N 
n 



Q 



m ^ u, u ^ m, K ^ n — K, q ^ q 
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Setting 

N^^, l^^^p' {=l-p), (5.24) 
we find, in the same way as formula (5.17) was gotten, that 



Pr{L = M) 



(1 -pOV""", (5.25) 



a different g- version of the classical Bernoulli distribution. 

6 g-contagious distribution 

Suppose we again, like in the preceding section, have an urn with m marked and n un- 
marked balls. We pick out one ball at random, record its value, and then return to the urn 
s + 1 balls identical to the one we just picked out. If s = 0, we return the ball itself, and 
this is the Bernoulli scheme; if s = — 1, we return nothing, and this is the hypergeometric 
scheme. For general s, the probability to pick out k marked (by '1') balls out of n draws 
is 



Pri^n = = (^) n2=o (rn + as) n^X' + / ^j=o {m + u + js). (6.1) 

This particular member of the family of the so-called "contagious distributions" was dis- 
covered by Eggenberger and Polya in 1923 [2, 8]. 
As a Qf-analog of this distribution, we set 



Pri^n = [«]) 
n 

K 



3.2) 



Obviously, for the case s = — 1 we recover the hypergeometric formula (5.2). For the case 
s = 0, we recover the classical Bernoulli formula (2.3) with p = [m]/[m + u]), not the 
g-Bernoulli formula (2.4). 

To justify formula (6.2), we need to check that 



K=0 



^ g(m+a«)(n-K)n^-i [^ + «5]n^-«-i[^ + ^s] = n!;-^[m + ^x + 7s]. (6.3) 



q-s 



To do that, we assume that s ^ and start with the formula (5.3) in the base Q = q 



E 

K=0 



'M' 




U 


K 


Q 


n — K 



Q 



(M-fc)(n-K) 



M + U 
n 



Q = 



(6.4) 
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which we rewrite as 



K=0 



Multiplying both sides by ([—fig)", setting 

M = —m/s, U = —u/s, 
and noticing that 

Q{M—K){n—K) q—s{—K—m/s){n—K) ^{m+Ks){n—K) 



(6.6) 



(6.7) 



we arrive at formula (6.3). The latter formula can be considered as a new g-analog of 
Newton's binomial. 

Similar to the hypergeometric case, we can arrive at the qf-contagious distribution (6.2) 
as a macroscopic object, 

6^ = Cl + ••• + Cn (6.8) 

from the microscopic formulae 

= qEt=oi"^+s^)ai^} + n^-^"^ [u + ps] / H^^Z^ [m + u + ^s]; (6.9) 



the latter formulae are suggested by the extreme cases k = n and /c = of formula (6.2): 

(6.10a) 
(6.10b) 



Pr{in = [n]) = n^Z ^^^Jy N = u + m, 



To verify that microscopic formulae (6.9) imply the macroscopic formula (6.2), we need 
to check that 



n 



\a{K)\=n—K 

Now, for the LHS of formula (6.11) we get: 



^ gEz"=o('n+«^Wi) = g'raCn-K) ^(-g-s^-E^^W [by (2.35)] 



(6.11) 



\a{K)\=n—K 

_ ^rn{n-K) i^^-s^-{n-K)K 



{m+SK){n— k) 



n 



(6.12) 



and this is exactly the RHS of formula (6.11). 
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7 q-uniform distribution 

A classical random variable X taking M + 1 discreet values vq < ... < vm, each with 
equal probability 1/(M + 1), represents a discrete uniform distribution. The values v'^^s are 
immaterial and can be taken as Vi = i, oi Vi = a + hi, oi Vi = [i], ... 
As a g-analog of this distribution, we set 

Pr{X = i)= q'/[M + 1], < i < M, (7.1) 

or 

Pr{X = [i]) = q'/[M + 1], < i < M, (7.2) 
so that 

<^>=«, (7.3) 

V^riX)^m^EL. ,,5, 

Consider n independent identically distributed, via the discrete uniform distribution, 
random variables Xi, The range of these variables is the quantity 

r„ = maa;(Xj) - mm(Xi), < T < M. (7.6) 

i i 

The random variable r„ has the following distribution ([4]) p. 240): 

Pr(r„ = 0) = 1/(M + l)"-\ (7.7a) 
Pr{rn = i) = {{£ + 1)'' - 2r + {£- 1)")(M + 1 - £)/(M + I)'', 1 < £ < M. 

(7.7b) 

This distribution is certainly different from those appearing in the preceding sections. 

Let us calculate the g-analog of the distribution (7.7). Taking as our basic definition 
formula (7.1), we have: 

M I 

Pr{rn = 0) = ^(Pr(X = i))" = [M + l]^n [M + 1]"; (7.8) 



=0 



^ M-l " / \ 



i=0 k=0 

M-l 



1 ^"^ [^]n«([2]"-[2]„n) 
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finally, for £ > 2, 

1 ^^"^ / \ 

J j=0 ^ all k's k(0)=0 K{e)=0 k(0)=k(^)=0 ^ 

1 ^^-^ 

= uTTm E ((9* + - + 9'''')" - (9'^' + - + I'^'y 

[Af + 1 - £] n / N 

For £ = 1, formula (7.10) reproduces formula (7.8). Thus, formulae 



[M + 1] 

Pr(r„ = 0) = ^, (7.11a) 

V J [^j + 1]" ' ^ ^ 

[M + 1 - £] n / X 

^^(^" = ^) = [jvf + I]'- + " t^^" + ^"^^ " J ' 1 ^ ^ ^ ^' 

(7.11b) 

are g-analogs of formulae (7.7). Notice, that 

Pr(ri > 0) = 0, Vg. (7.12) 

Although not immediately apparent, formulae (7.11) are not the only g-analogs of the 
classical formulae (7.7). For example, for n = 2, formulae 

= "> = ImTII' 



are different from formulae (7.11) 



n=2 



8 Concluding remarks 

The approach to g-Probability taken in this paper leaves the rules of classical probabil- 
ity intact and only (j'-dcforms some basic probability distributions. It is quite likely that 
one can develop some new/bizarre rules of g-probability which contradict the comfort- 
ably familiar intuition, similar to what Quantum-mechanical interpretations appear to 
a Classical-mechanical disciple. The most direct route to such new rules probably goes 
through basic continuous probability distributions when integral J {■)dx is replaced by the 
g- integral j {■)dqX. 
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